arXiv:1502.02082v2 [quant-ph] 12 May 2015 


Maximally genuine multipartite entangled mixed 
X-states of A^-qubits 


Paulo E M F Mendonga^ Seyed Mohammad Hashemi 
RafsanjanP, Diogenes GalettP and Marcelo A MarchiollP 

1 Academia da Forga Aerea, C.P. 970, 13.643-970 Pirassununga, SP, Brazil 

2 Rochester Theory Center and the Department of Physics & Astronomy, University 
of Rochester, Rochester, NY 14627, USA 

3 Institute de Fisica Teorica, Universidade Estadual Paulista, Rua Dr. Bento 
Teobaldo Ferraz 271, Bloco II, Barra Funda, 01140-070 Sao Paulo, SP, Brazil 

4 Avenida General Osorio 414, centre, 14.870-100 Jaboticabal, SP, Brazil 

E-mail: pmendonca@gmail.com, hashemi@pas.rochester.edu, 
galetti@ift.unesp.br, marcelo_march@bol.com.br 

Abstract. For every possible spectrum of 2^-dimensional density operators, we 
construct an 7V-qubit X-state of same spectrum and maximal genuine multipartite 
(GM-) concurrence, hence characterizing a global unitary transformation that — 
constrained to output X-states — maximizes the GM-concurrence of an arbitrary input 
mixed state of A-qubits. We also apply semidefinite programming methods to obtain 
A-qubit X-states with maximal GM-concurrence for a given purity and to provide an 
alternative proof of optimality of a recently proposed set of density matrices for the 
purpose, the so-called X-MEMS. Furthermore, we introduce a numerical strategy to 
tailor a quantum operation that converts between any two given density matrices using 
a relatively small number of Kraus operators. We apply our strategy to design short 
operator-sum representations for the transformation between any given A-qubit mixed 
state and a corresponding X-MEMS of same purity. 
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1. Introduction 

In the framework of quantum information theory, mixed A-qubit X-states synthesize 
a family of quantum states whose inherent correlations are much easier to quantify 
than is generally the case. The prefix ‘X’ is motivated by the shape of their density 
matrices written in the computational basis [^, whose non-zero entries are either 
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diagonal or anti-diagonal (or, otherwise, can be brought to this form via a local unitary 
(LU) transformation). Owing to this sparse structure, that includes important states 
(e.g.. Bell’s (^, Werner’s [^, isotropic [^, GHZ [^, etc.), analytical investigations of 
entanglement properties (^|6 12 and quantum discord [l3}|^ in Wqubit X-states have 
lately become an active and fruitful field of research. 

Retrospectively, important members of the class of two-qubit X-states were 
identified in 26,27 , where the concept of maximally entangled mixed states (MEMS) 
was introduced and characterized. From these seminal works, worthy of note is the 
observation that two-qubit states of a fixed spectrum and maximal entanglement (as 
measured by concurrence, negativity or relative entropy of entanglement) can always be 
found in the X-form. Subsequently, Munro et al. characterized two-qubit states 

of maximal entanglement for a fixed mixedness (as measured by purity, linear entropy 
or von Neumman entropy), once again obtaining X-states as results. 

In spite of these early achievements, to date, little has been accomplished in 
extending the characterization of MEMS beyond two-qubits. Largely, this is because 
sensible measures of genuine multipartite entanglement have been identified only 
recently 30,31 and are generally hard to evaluate, let alone maximize. 

A first important step toward the identification of iV-qubit MEMS for N > 2 
was given by Hashemi Rafsanjani et al. [^, who showed that the GM-concurrence of N- 
qubit X-states admits a simple closed formula, amenable to maximization. Although the 
resulting optimal states of this maximization cannot be guaranteed to be actual MEMS, 
at least they are provably MEMS amongst all Wqubit X-states. Therefore, in 


10 


Agarwal and Hashemi Rafsanjani maximized the X-state GM-concurrence formula under 
the constraint of a fixed linear entropy, determining the so-called X-MEMS. 

In this paper, we enlarge the scope of the term X-MEMS to enclose two classes 
of X-states: X-MEMS with respect to (wrt) spectrum, referring to those Wqubit X- 
states of maximal GM-concurrence for a fixed spectrum, in analogy to the original 
MEMS introduced in 26,27 ; and X-MEMS wrt purity, referring to Wqubit X-states 


29 


that maximize the GM-concurrence for a fixed value of purity, in parallel with 
Our main results initially consist of: (i) a complete characterization of X-MEMS wrt 
spectrum, and (ii) a demonstration that X-MEMS wrt purity can be obtained from 
the solution of a semidefinite program (SDP) [^,33 , by which means (iii) we provide 
an alternative proof of optimality of the states obtained in 10 . Moreover, 


IV we 


characterize the unitary transformation that maximizes the X-state GM-concurrence 
formula of an arbitrary Wqubit state, generalizing the result of for N = 2. Finally, 
of independent interest (but also relevant in this context), (v) we construct a family 
of iterated SDRs whose solutions produce quantum operations (GPTP maps) that 
implement a desired state transformation with a decreasing number of Kraus operators. 
The method is illustrated with the determination of short operator-sum representations 
for the conversion between an arbitrary input state of purity P and a corresponding 
X-MEMS wrt purity P. 

Our paper is structured as follows. In section]^ we briefly review the concept of GM- 
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concurrence and, in particular, its simple formula for iV-qubit X-states. In section we 
characterize X-MEMS wrt spectrum and the unitary transformations that produce such 
states from arbitrary A^-qubit density matrices. In section X-MEMS wrt purity are 
constructed and have their optimality established via SDP theory. Section [^outlines the 
iterated SDP method to design quantum state transformation with few Kraus operators 
and exemplihes the method while producing X-MEMS wrt purity from arbitrary input 
states of the same purity. Finally, section summarizes our results and discusses some 
possible avenues of future work. 


2. GM-concurrence of A^-qubit X-states 


In this section, we present the formula for the GM-concurrence of X-qubit X-states 
obtained in 1^. For the beneht of the reader unfamiliar with the current literature on 


multipartite entanglement (in particular, [9,31,^), we start by reviewing some key 
dehnitions concerning A^-qubit X-states, GM-entanglement and GM-concurrence. 

To begin with, let us introduce some notation. Throughout, denotes the 
(complex) Hilbert space of dimension di, whereas i3(HdJ denotes the set of (bounded) 
linear operators acting on H^.. The set of all possible bipartitions of {1,2,... ,N} is 
denoted by P and a particular bipartition in P is denoted by P^ (with rj ranging 

from 1 to 2^“^ — 1). Partial traces over Hilbert spaces whose labels i belong to 
are concisely indicated as trn . 


Definition 1 An operator px € H(H 2 ® ® H 2 ) represents an X-qubit X-state if and 

only if, in the computational basis {|bini)}^^Q“^ (and up to a LU-transformation), it 
assumes the matrix form 

ai 

02 r 2 e *'^2 


an 

rp p Vs 


( 1 ) 


rie 


-i0l 


r^e *'^2 


&2 


bi 


where n := 2^ ^ and, for every integer k G [l,n], we have Ok, hk, Vk G M+, (pk G [0,27r], 

n 

5^(afc + bk) = 1 and 0 < Tfc < y /Ukbk ■ (2) 

k=l 

While ([^ visually justihes the prehx ‘X’, the conditions (|^ ensure the normalization 
and positive semidehniteness of px- As a glance at ([^ demonstrates, the index k G [1, n] 
can be regarded as a label for uncoupled bidimensional subspaces. That any X-qubit 
X-state is decomposable into n such subspaces is a key property that will be implicitly 
exploited throughout this paper. 
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Although we are only interested in the entanglement properties of A^-qubit X-states, 
we proceed with a general definition of GM-entanglement. 

Definition 2 An N-partite density operator p G ® ® ® is GM- 

entangled if and only if it is not biseparable. 


To understand the concept of biseparability, consider first its definition for pure states. 

Definition 3 An N-partite state 1-0) G 0 ® ... 0 is biseparable if and only 

if there is a Hilbert space bipartition (g) 0 0 ... 0 and a pair of 

states IV’a) G I^Jb) ^ such that \^p) = \^pA) ® IV’s)- 

Note that definition implies that a biseparable state is not necessarily separable, as 
there might be entanglement within and/or H^. It then follows from definition]^ 
that the condition for GM-entanglement is generally more stringent than the condition 
for bipartite entanglement, for example. In fact, GM-entanglement only occurs when 
bipartite entanglement is observed across all possible bipartitions of 

The notion of biseparability is extended to mixed states as follows. 

Definition 4 An N-partite density operator p G 0.. is biseparable 

if and only if it can be decomposed in an ensemble of biseparable pure states, that is 

p = ^Pi\A){A\, (3) 

i 

where YhiPi — ^ iV’i) biseparable (even if with respect to different bipartitions 

of Hdi 0 0 ... G) 

The above definitions provide a formal criterion to determine whether a general 
mixed state is GM-entangled or not. A further step was given by Ma et al. [^, who 
introduced the GM-entanglement measure named GM-concurrence. 


Definition 5 The GM-concurrence of an N-partite pure state \^fl) G . -(gHrfjY 

is given by 

Ccrnim ■■= mm_^_^^ - tr[pij , (4) 

where pa,, '■= trB,,[|V’)(V’|]- For N-partite density operators p G 
the GM-concurrence is obtained via the convex roof construction 

C'gm(p) = inf VpiGGM(IV’i)), (5) 

{pi.m)} 

I 

with the infimum taken over all possible ensembles {pi, Ifi)} that realize p. 


The GM-concurrence takes its name from the fact that, in the case of two-qubit systems, 
it matches the Wootters concurrence 35 and, more generally, can be shown 31 to 


satisfy the following minimal requirements for any GM-entanglement measure: 
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GM-entanglement detection: 

Cgm(p) > 0, 

with saturation if and only if p is biseparable. 
Convexity: 


( 6 ) 


G 


GM 


PiCGMiPi) ■ 

\ i / i 


(7) 


( 8 ) 


Monotonicity under local operations and classical communication (Clocg): 

Cgm(Clocg[p]) < C'gm(p) • 

Invariance under LU-transformations (17 l): 

Cgm{UlPUI) = Cgm{p) ■ (9) 

Though well motivated, Cgm(p) is generally hard to evaluate due to the infimum 


over all ensembles that realize p. To alleviate this problem, the authors of 31 relied on 


certain sufficient criteria for GM-entanglement detection proposed by Huber et al. 34 


to determine computable lower bounds for Ggm- In particular, if the main- and anti¬ 
diagonal entries of p are parametrized as in Q (the remaining entries being arbitrary), 
then one of Ma’s lower bounds reads (see Appendix A] for an explicit derivation) 


Cgm(p) > 2 max 0, max 

fce[l,n] 


( 10 ) 


Remarkably, as shown by Hashemi Rafsanjani et al. [^, this lower bound is saturated 
when p = px, namely. 


Cgm{Px) = 2 max 


0, max 


rk 


( 11 ) 


The fact that A^-qubit X-states have their GM-concurrence expressed as a closed 
formula cannot be overstated. It contrasts with the great difficulty involved in merely 
detecting GM-entanglement in more general systems, not to mention quantifying it. Of 
course, this result becomes even more appealing when one notices that A^-qubit states of 
practical interest do occur in the X-form (see, e.g., (^), or otherwise can usually be well 
approximated to it via LU-transformations [^. Finally, it is interesting that for GHZ- 
diagonal states (X-states with a* = fej), the value of GM-concurrence is proportional to 
the distance of the GHZ-state to the set of biseparable states (Ml. 


3. X-MEMS with respect to spectrum 

As mentioned before, the two-qubit MEMS with a given spectrum, characterized 
in (^j^, are X-states. In this section, we assume that this is also true in the A^-qubit 
case {N > 2), and characterize the ‘W-qubit MEMS” resulting from this assumption. 
Since it is not known in which circumstances the restriction to the set of X-states is an 
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active constraintj|] for N > 2, we adopt the nomenclature introduced in 10 and talk 
about X-MEMS instead of simply MEMS. 

The results of this section are summarized in Theorem which is deliberately 


presented in close resemblance to the statement of the related Theorem presented in 27 


regarding the case N = 2. The proofs, however, are established in very different ways. 

Theorem 1 The maximal GM-concurrence attainable by an N-gubit X-state of 
spectrum A, determined by the eigenvalues Xi > X 2 > > X 2 n (with n := 2^~^), 

is given by 


max 


0, Ai — An+i — 2 a /XeX 


2n+2-£ 


e=2 


( 12 ) 


Moreover, any N-gubit density matrix p\ (of spectrum A ) can be coherently transformed 
into p'f^, an isospectral N-qubit X-density matrix of maximal GM-concurrence, according 
to = lA.pifil\ with the unitary U given by 

N 

U=\ 5^ [/, 1 . (13) 



-1 

1- 

(M 

_1 

V22 


In (13), the following definitions apply: {Uk}^=i is a set of arbitrary single qubit unitary 
operations, is an arbitrary unitary and diagonal matrix, $ is the unitary matrix 
formed from the eigenvectors of p\ (such that p\ = and 

n ^ ^ n—1 

Vn = V 12 + Ei,i , V 12 = , V 21 = —^Enp , V 22 = — V 21 + • (14) 

i=2 v2 V2 

Here, Eij is the n-dimensional matrix whose only non-zero entry is equal to 1 and 
occupies the ith row and jth column. 


An immediate remark is that, as expected, both the optimal GM-concurrence (12) 


and the unitary transformation (13) reduce to the corresponding expressions in 27 


Theorem 1] when N = 2. The remainder of this section is devoted to proving the 
theorem for N > 2. Essentially, our proof consists of a direct maximization of the 
GM-concurrence formula of A-qubit X-states under the constraint of a fixed spectrum. 

Let denote a generic A-qubit X-density matrix of spectrum A. We start by 
taking matrix (j^ as a parametrization for and writing a general formula for Ggm(p(;^) 
in terms of its eigenvalues 




Ofc -|- bk 


± 


+ d\ for every k G [l,n] 


(15) 


Here, A^ denote the greatest (-I-) and smallest (—) eigenvalues of p'^ associated with 
the bidimensional subspace labelled by k, and dk ■= {bk — ak)/2. It follows trivially from 


(15) that. 


+ dl = 


Xt -XI 


and 




( 16 ) 


§ By an active constraint we mean a restriction that is not satisfied unless it is explicitly imposed. 
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which used in ( p!I| yields 

Cgm(Pa) = 2 max 



A+-Ar 




, + 


i=2 


( 17 ) 


where, without loss of generality, we fixed the label /c = 1 to the subspace whose value 
of rk - j is the largest. 

Our goal now is to maximize 0 under the contraint that the set of eigenvalues 
{^f}'j=i matches the set of eigenvalues of the arbitrary (but given) iV-qubit density 
matrix pa, he., A = A. To do that, let us first consider the optimization over di and 
{^j}]= 2 - Although these variables are constrained by (j^ and related to by (15), 

we will momentarily ignore these contraints. By doing so, we significantly simplify the 
optimization procedure at the expense of risking over-maximization Ogm- Nevertheless, 
as we will soon demonstrate, the resulting maximal is actually attainable, meaning that 
our simplifying assumptions are harmless. With that in mind we set, for every j G [2, n], 


rj = di = 0, (18) 

which clearly maximizes 0 over Tj and di. 

At this point, we are left with the maximization over written as 

maximize Xf — u ■ v subject to {A^}Ei = {^j}‘j=i ! (19) 

where Ai > ... > X 2 n are the eigenvalues of the arbitrary (but given) iV-qubit density 
matrix Pa and the vectors u,v E are given by 



Here, we aim to assign to each variable in an eigenvalue of p\, in such a way 

that Aj*" is maximal and u ■ v is minimal. To maximize Afy we simply assign to it the 
largest eigenvalue of pa, i.e.. 


AJ = Ai. 


( 22 ) 


To minimize u ■ v, first notice that u and v display the same entries in the reversed 
order, with fyE" occupying the central position in both vectors. It follows from the 
rearrangement inequality (see, e.g., 39, Theorem 368, page 261]) that the scalar product 
between two vectors defined up to the ordering of their entries is minimized if and only if 
they are sorted in opposite directions. So, we make the entries of u and v monotonically 
increasing and decreasing, respectively, by assigning, for every j G [2,?7,], 

-^1 = A„+i, and A • = X 2 n+ 2 -j ■ (23) 


A+ = A,, 


Substituting the identities (18), (22) and (23) in (16) and solving the resulting 
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every j e [2,n], 


, we 


obtain thatjj]] for 


ri = 


Ai - A 


n+l 


ai = bi = 


Ai-A 


n+l 


Tj — 0 , Oj — Xj 


and bj = A 


2n+2-j ■ (24) 

IS 


2 ' 2 
By plugging (24) into 0 . it is easily seen that ([T^ holds. In order to see that ( [l^ i 
physically attainable, substitute (24) into matrix ([^ (and set = 0 for every k G [1, n]), 
to get 


Pa = 


Ai + A„+i 


Ai-A 


n+l 


2A, 


0 


2A„ 0 

0 2 A „+2 


0 2A2n 

Ai — A,i+i Ai + A,i+i 

It is immediate to check that is a valid X-density matrix with the same spectrum of 


(25) 


Pa and GM-concurrence given by (12). 


Finally, let us establish (13). Since p\ and p'^ are isospectral, we can write 

Pa = $A$'^ and p'^ = VKV^, (26) 

where $ is the matrix of the eigenvectors of Pa and V is the matrix of the eigenvectors of 


P(y given by (25). Some simple linear algebra shows that, for A = diag[Ai, A 2 ,..., A 2 n], 


the matrix V admits the block decomposition specihed in (14). Thus, combining the 


two identities in (26) to eliminate A, we arrive at 


p(v = UptfU ), where tl = 


Vn 

V12 

V21 

V22 




(27) 


As noted before, the X-MEMS of (25) are unique up to LU-transformations, for 


which reason, in (13), the expression of lA. appears pre-multiplied by an arbitrary LU- 


transformation. Furthermore, for sake of generality, we have also multiplied an arbitrary 
(generally non-local) diagonal unitary matrix T)^ in (13). It should be emphasized, 
though, that has obviously no effect on the output state p'^. 

Let us conclude this section by answering a central question that arises from the 
present work: are iV-qubit X-MEMS actual A-qubit MEMS? Although this is long 
known to be the case for N = 2 27 , indications that the same may also hold for 
N = 3 have only recently appeared in the work of Hedemann 40 . Alas, to the best 


of our knowledge, the topic seems to be utterly unexplored for A > 4. To see that 
II As a matter of fact, many other solutions can be obtained by interchanging the values of Uj and bj 


indicated in (24) for any j S [2,n]. However, this does not lead to essentially new X-MEMS, since the 
X-MEMS corresponding to these solutions can always be generated from the X-MEMS corresponding 


to (24) via a LU-transformation. 
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that conjecture cannot hold in general, note that an affirmative answer (combined with 


(12)), would imply that A^-qubit density matrices whose eigenvalues Ai > A 2 > ... > X 2 n 


satisfy 


n 

Ai < An+i + 2 \ \/XiX 


e=2 


2n+2-e 


(28) 


cannot acquire GM-entanglement by means of a global unitary transformation. However, 
as recently shown by Huber et al. [^, A^-qubit thermal states of arbitrarily high 
temperatures and N sufficiently large (represented, in the computational basis, by 


diagonal density matrices arbitrarily close to the identity, hence fulhlling (28)) can 


acquire GM-entanglement by means of rotations to Dicke-like (non-X) states, thus 
providing a counter-example to the original conjecture. 


4. X-MEMS with respect to purity 

Since the purity of a given quantum state can be expressed as the sum of the squares 
of their eigenvalues, there are, in general, many spectra that realize a hxed value of 
purity P. In this section, we consider the problem of searching for a maximizer of the 
X-qubit X-state GM-concurrence formula amongst all A^-qubit spectra that realize P. 
Any X-state with such an optimal spectrum is referred to as an X-MEMS wrt purity, 
denoted by p'p. 

As a matter of fact, X-qubit X-MEMS wrt purity have been recently determined in 


the work of Agarwal and Hashemi Rafsanjani 10 , whose main Endings are summarized 
in the statement of the following theorem. 

Theorem 2 (Agarwal and Hashemi Rafsanjani) For any value of purity P G 
]l/(n-|- 1),1] with n := , the maximal GM-concurrence attainable by an N-gubit 

X-state of purity P is 27 , where the parameter j g]0, 1/2] is determined by P according 
to 


7 : = 


1 


2(n+l) 






n+l 

71+3 

(n+l) 


n+3 

(n+l)^ 

2<P<1 


(29) 


Up to LU-transformations, every N-qubit X-density matrix of purity P that achieves 


maximal GM-concurrence is given by (25) with 

An+l = /(7) 


7 


and 


'^n+j 0 


(30) 


-''i = /(7)+7. -'j=s(7). 

for every j G (2,n], being f and g defined as follows: 

7 ^ ^f 0 < 7 < 7 ^ 

7 if I 

In what follows, we give an alternative proof of this theorem by exploiting the 
results of Theorem and the theory of semidefinite programming 32,33 . 


fix) ■= 


and g{'j) : = 


1 - 2/(7) 


n 


(31) 
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Since purity is determined by the spectrum (and not the other way around), the 
specihcation of a spectrum generally represents a stronger constraint than a specihcation 
of purity. As a result, if p'p is a X-MEMS wrt purity P and has spectrum Ap, then p'p is 
also a X-MEMS wrt to the spectrum Ap. In other words, every X-MEMS wrt purity can 
be regarded as a X-MEMS wrt some spectrum, in which case every p'p is of the form 


(25) up to a LU-transformation. Thanks to this, we can determine p'p by maximizing 


the GM-concurrence of (25) over all sets of physical eigenvalues that realize P, which 


yields the optimization problem]^ 
maximize Ai 


subject to A 


An+l ~ 


> As > 


A/AjA2n+2-i 

i=2 


2n 

> Asn > 0 , Afc = 1 and A 


2n 


1 ^ 


k=l 


2 _ p 

fc ~ • 


(32) 


k=l 


Next, we give a few arguments that allow some simplification of problem (32). 
First, there is no need to explicitly require the ordering Ai > As > ... > Asn since 
we know in advance (rearrangement inequality) that this particular ordering is a 
necessary condition for the maximization of the considered objective function and will 
be thus satisfied anyway. Second, the objective function is clearly maximized if we set 
Asn+ 2 -j = 0 for every j G [2,n], which does not violate any problem constraint and 
reduces the equality constraints to 


An+i 1 ^ ^ A/j and | 1 ^ ^ A/j | -|- ^ ^ A 


2 _ p 


(33) 


k=l 


k=l 


k=l 


Given these two points, we can replace the original (non-linear) objective function with 
the linear function Ai — A„+i, and all the inequality constraints with a single one: 
An+i P 0. Finally, without altering the solution of the problem, we can replace the 
equality in the quadratic constraint with the inequality ‘<’, thus enlarging the set of 


feasible points to its convex hull 42, Ghapter 32]. Accordingly, we end up with the 
equivalent optimization problem on n real variables: 


maximize 


1 -|- 2Ai -|- ^ ^ Aa 


fc =2 


subject to Afc < 1 and ( 1 — 


(34) 


k=l 


k=l 


k=l 


As explained in Appendix A, the quadratic (convex) constraint in (34) can be 


written as a linear matrix inequality (LMI), turning problem (34) into the SDP (A.7). 


^ Note that here {Xk}k=i represents the set of optimization variables, as opposed to a fixed set of 
eigenvalues (as it was the case in Sec.j^. 

+ Up to (irrelevant) permutations of the form Xj O Asn+s-j for any j € [2,n]. 
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It is straightforward to see that ( |A.7[ ) admits the SDP standard ineguality /or? 7 i[^ 

n 

minimize c'''A subject to Fq + FiXi > 0 (35) 

with 


2=1 



Qn 




On 



c :=-b, Po := 


P- 1 


and Fi := 


2 




1 



-1 


(36) 


where A, b, Qn were defined in the Appendix A and G denotes the n- 


dimensional vector whose only non-zero entry is 1 and occupies the Ah row. 

The fact that the design of the maximally GM-entangled X-states of A^-qubits 
can be cast as a SDP is interesting on its own. First, SDPs are convex optimization 


problems 33 and, as such, have the desirable property that any local optimum is 
necessarilly a global optimum. Second, many efficient numerical methods have been 


devised to solve SDPs (for a review, see, e.g., 32 and references therein). These methods 


have excellent convergence properties and output a “certificate of convergence”, i.e., an 
interval within which the optimal value of the objective function must lie. Typically, 
with no more than 30 iterations, this interval can be made arbitrarily small. Third (and 
most importantly for our purposes), a powerful duality theory exists for SDPs and can 
be employed to rigorously prove the optimality of an ansatz solution. 

Before proceeding with our proof, let us exploit the aforementioned numerical 
virtues of SDPs to provide a first evidence of the optimality of the GM-concurrence 27 
(cf. (29)) and of the spectrum (30). In figure we plot these analytical expressions 
(lines) along with the numerical solutions (symbols) of problem (35), obtained by 


running the MATLAB-based solver SeDuMi for several combinations of P and N. 
In our numerical computations, we set SeDuMi’s precision to 10“^^, which establishes 
the largest acceptable length of the aforementioned “error interval”. 

Plot shows the converged numerical values of the objective function along with 
an analytic plot of 2'j{P). The difference between the numerical and analytical values is 
found to be of the order of 10“^^, which strongly suggests that 27 (P) is the maximal N- 
qubit X-state GM-concurrence wrt purity. Plots in turn, indicate the optimality 
of (30). In particular, plot & reveals an excellent agreement between the converged 
numerical values of Ai and the function /( 7 (P)) -|- 7 (P). Similarly, plot illustrates 
the coincidence between the converged numerical values of X 2 ,...,n (which resulted all the 


same up to numerical precision) and the function g{'y{P)), cf. (30) and (31). 

We now briefly review an important result from the SDP duality theory that will 
be subsequently used to establish the optimality of (30). For a thorough account on 


It should be noted that, although problems (A.7) and (|35|) are solved by the same set of eigenvalues 


the resulting optimal values of the two problems are not exactly the same. That is because, to arrive at 


problem (351, we removed the summand —1 from the objective function of (A.71 (and used the property 


max^, a(a;) = — min 2 ;(—a(a;)) to replace the maximization with the minimization). As a result, if tt* 


and p* denote the optimal value of problems (A.7) and (35), respectively, then tt* = —1 —p 
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a 





Figure 1. Matching between the numerical solution (symbols) of the SDP (35) and 
the corresponding analytical formulas (lines) for the optimal GM-concurrence and 
spectrum (cf. theorem|^. Each plot considers the parameter values iV = 2,3,4, 5 and 
a uniform sampling of P G ]l/(n + l), 1] . The legend in plot (a) also applies to plots (b) 
and (c). The observed agreement between the lines and symbols strongly suggests the 
optimality of the analytical formulas in the statement of theoremFrom left to right, 
the unlabelled tick marks correspond to the purity value := (n + 3)/(n + 1)^ for 
n = 16,8,4,2 (or, equivalently, N = 5,4,3, 2). Remarkably, although all the plotted 
functions are continuous, they fail to be smooth: the second derivative of Cgm(P) and 
the first derivatives of are discontinuous at P = Pn- 


this theory, we refer the reader to Chapter 5]. The dual problem of the SDP (35) 
(henceforth called primal problem), is another SDP given by 

maximize — tr[i^o-^] subject to Z > 0 and tr[i^jZ] = Cj \/i = 1,... ,n . (37) 

Here, the variable to be optimized is the matrix Z, whereas the vector c and the matrices 
-Fq and Fi are the same as the ones appearing in the primal problem (in our particular 
case, dehned in (36)). Let p denote any feasible value of the primal problem (35), and 
denote by p* its optimal value. Similarly, let d and d* denote, respectively, any feasible 
value and the optimal value of the dual problem (37). It is obvious that p > p* and 
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that d* > d. Less obvious — but also tru^— is that p > d ioi every primal and dual 
feasible values (in particular, p* > d*)/m. such a way that 

P > P* > d* > d, (38) 

a general result known as weak duality. Next, we show how these inequalities can be 
used to prove the optimality of our ansatz solution. 

As some straightforward computation shows, the spectrum (30) can be initially 
regarded as a primal feasible point that yields a primal feasible value p = —1 — 2'y. 
However, if we can hud a particular dual feasible point that yields a dual feasible value 
d = p, then weak duality implies that p = p*, meaning that (30) is, indeed, a primal 
optimal point. We claim that the block matrix Z, described below, provides such dual 
feasible point: 


S = 


3i 

V3i32j1-i 

-1-33 + 1 


^3132^71—1 

32+-1 

(-| - 33 + 7 ) 3n-l 


-1-33 + 1 

(-| - 33 + t) jn-l 

33 



34 


(39) 


where jn G 


Dnx 1 


and Jn G 


denote the “all-one” n-dimensional vector and matrix. 


respectively. In addition, 

if 


3i : = 


2 + 27 + ^ 

(l-n)(l+27)2 

2 {l- 2 n 7 ) 

(1+7)" 


0<7<^ 


if 

71+1 


<7<i 


32 : — 


27 

(n- 27 )" 

2(1— 7 l)(l—2717) 


if 

if 


0<7<;^ 


n+l 


< 7 < t 


33 


34 := 


27 

l—n 

2 (l- 2 n 7 ) 


if 0<7<^^ 
if 7TT<7<^ 


(40) 


1 + W 


1-27 


2717 


if 0 < 7 < 

if WI<7<f 


To verify our claim, we hrst show that Z establishes d = p. Indeed, d = —ti [FqZ] , 


which can be evaluated with the aid of (36) and (39) to give 

1 


d = 


n + l 


[ 3 in + 2 (n - 1)(32 - v^i 32 )] - h{P - 1 ) - 34 • 


(41) 


Then, plugging (40) into (41) and expressing P as the sum of the squares of the 
eigenvalues in (30), we obtain d = —1 — 27 = p, as claimed. 

Now, we show that Z satishes the constraints of problem (37). Regardless of the 
values 3 i, 2 , 3 , 4 , it is easy to check that tr [FiZ] = —2 and tr [i^ 2 ,...,n- 2 ] = — 1 , as required. 
Finally, the condition Z > 0 can be checked by an explicit study of the eigenvalues of 
Z. Simple inspection of (40) shows that 34 > 0. Furthermore, with some cumbersome 

ji To see that, note that the constraints of the primal and dual problems allow us to write: 


I — d = c^X + tr [FqZ] = ^ tr [F^Z] + tr [FqZ] = tr 




> 0 . 
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simplification procedure, the characteristic polynomial of the first (n + l)-dimensional 
block of 2 takes the form 


X 


- An(7)0:" = 0 


(42) 


with 


An(7):= 


(n+l)(l+ 27 )+(n+ 3 ) 7 ^ 

27 

n^+ 2 (n— 1 )— 47 + 4 n 72 


if 0 < 7 < ^ 

if 


(43) 


2 (-l+ 2 n 7 ) 

thus establishing x = An( 7 ) as the only non-zero eigenvalue of this block. Clearly, the 
first branch of (43) is strictly positive. To see that this is also true for the second branch, 
define 


Tn( 7 ) := 2(?t. - 1) - 47 4n7^ and Gnij) ■= 2{-l2n'y) (44) 


in such a way that, for 7 G [l/(u + 1), 1/2], An(7) = Fn{'j)/Gn{'j)- We conclude our 
proof by showing that, for 7 > l/(n -f- 1), both Fni^j) and Gn{'y) are strictly positive, 
hence so is An( 7 ). To establish the positivity of Gn(7), note the following implications: 


7 > 


u -|- 1 


7 > ^ ^ G„(7) > 0. 


(45) 


The first implication follows from the fact that, for the relevant values of n (recall that 
u > 2), the inequality n -|- 1 < 2n holds trivially. The positivity of Fn{'j), in turn, can 
be established by noting that 




n -I- 1 


(n — 1)(?T, -I- 3)[2 -I- n{n -|- 2)] 
{n 1)2 


> 0 and 


d'y 


= 2G„(7)>0. (46) 


So, we see that Fni^j) is already positive at 7 = l/(n -|- 1) and monotonically increasing 
for 7 > l/(n -I- 1 ). 


5. Converting between density matrices with few Kraus operators 

For any given pair of isospectral density matrices, it is possible to find a unitary 
transformation that maps one density matrix into the other. In Theorem for 
example, all unitary transformations that map an arbitrary A^-qubit density matrix 
into a corresponding X-MEMS wrt spectrum were explicitly constructed. However, if 
two density matrices have different spectra (e.g., an arbitrary density matrix and a 
corresponding X-MEMS wrt purity, cf. Sec. |^, then there is no unitary map capable 
of converting between the two, in which case one must resort to more general quantum 
operations to implement the desired state transformation. In this section, we introduce 
a numerical scheme to design a quantum channel C, as modelled by a completely positive 
and trace preserving (CPTP) map, that promotes the conversion between any two 
given density matrices. Furthermore, our scheme constrains the resulting map to be 
“economical” wrt certain resources utilized in its implementation. 
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To understand what economical means in this context, consider the general 
representation of C in terms of Kraus operators G -B(Hd): 

Cip) = yM^pMl with yMlM^ = l. (47) 


m 


m 


A trivial choice of C, such that p = C{p), is the channel that maps every density matrix 
to ~p, which admits the following minimal set of Kraus operators: 

= .(48) 

where d := dimp, r := rankp, {|z^)}y=i is an orthonormal basis formed from the 
eigenvectors of p and are the non-zero eigenvalues of p (corresponding to the 

eigenvectors |/i)), i.e., 


P = 




(49) 


/ 4=1 


Practically, though, implementing (48) can be considered an overkill. In fact, since we 


only require p h->■ p, it might be possible to find a quantum channel that implements the 
desired state transformation with a number of Kraus operators much smaller than dx. In 
other words, there might be a more economical CPTP map for the state transformation 

p^p. 

With that mind-set, consider the task of determining, amongst every CPTP map C 
that satisfies C(p) = p, those that can be decomposed with the smallest possible number 
of Kraus operators. Mathematically, this leads to an optimization problem that can be 
nicely expressed with the aid of the Choi-Jamiolkoswki isomorphism 44 47 , which 
brings CPTP maps C ; BiH^) —> BiH^) into a one-to-one correspondence with positive 
semidefinite matrices C G i3(Hd® Hrf) such that tr 2 [C!:] = Id (here and throughout, tr^, 
denotes the partial trace over the xth d-dimensional subsystem). In particular, given a 
CPTP map C : B{V\d) —)■ i3(Hrf), its corresponding Choi-Jamiolkoswki matrix is 

(!:= (X®C)|T)(T| , (50) 

where X is the identity map on B{V\d) and |\k) is the (unnormalized) maximally entangled 
state |J/) = ® {|h^)}a=i a fixed orthonormal basis for H^. In this 

framework, the minimal Kraus decompositions of C can be shown to have rank C Kraus 
, and the equation C(p) = p is equivalent to tri [(p"*” (g) 1^)6!:] = p (the 

Accordingly, the 


operators 


48 


transposition being taken wrt p written in the basis {|h^)}„^;^) 49 
optimization problem takes the form 


minimize 


subject to 


rank C 

C > 0, 


tr 2 [e:] ^ Id and tr, [(/a^ ® ld)e:] ^ p. 


(51) 


Problem (51) is an example of rank minimization problem (RMP) with SDP 


constraints. Although special cases of this problem have been solved (see, e.g., and 
references therein), RMP are, in general, computationally intractable (NP-hard) due 
to the non-smoothness and non-convexity of the rank function. Fortunately, though. 
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heuristic methods exist to efficiently approximate their solutions 51 53 . Essentially, 


these heuristics rely on replacing the rank with a surrogate function, in such a way that 
the resulting problem can be handled with standard SDP solvers. Next, we consider the 


application of two such methods to problem (51). 


5.1. Trace heuristic 

Basically, it consists of replacing rankC with tr[C]. Intuitively, the replacement makes 
sense because sparse vectors tend to have a £i-norm smaller than dense vectors. So, 
by forming a vector from the eigenvalues of C and minimizing its f'l-norm (which 
corresponds to minimize tr[C] when (t > 0), we may be effectively vanishing some 
eigenvalues of C and, hence, reducing its rank (see 51,52 for a more technical 
justihcation of the trace heuristic in terms of the convex envelope of the rank). To 
a large extent, the appeal of the trace heuristic comes from the fact that it provides a 
linear objective function for the optimization problem, which usually means that the 
solution can be efficiently obtained (at least numerically). Unfortunately, in the case of 


problem (51), the constraint tr 2 [C] = 1^ implies that tr[C] = d, so no minimization can 
actually occur. As a result, the trace heuristic is useless for our purposes. 


5.2. Log-det heuristic 


Introduced in 51,53, the log-det heuristic can be considered a rehnement of the 


trace heuristic. Consists of replacing lank^ with logdet(C!: -|- dl^a), where d > 0 is 
a regularization constant. The value of S can be made arbitrarily small, and it is used 
to avoid an ill-dehned objective function as C gets singular (det C = 0). Intuitively, it is 
expected that, for C > 0, the determinant det(C-|-dlrf 2 ) will decrease as the eigenvalues 
of C vanish. So, we attempt to minimize the function log det (C-1-51^2), which plays the 


role of a smooth and concave 33 surrogate of the rank. 


The optimization problem resulting from the application of the log-det heuristic 


to (51) is a minimization of a concave function over a convex set, and it is thus non- 


convex. In order to obtain a related convex optimization problem, the log-det objective 
function is expanded (to the hrst order) in a Taylor series about a hxed d^-dimensional 
matrix Cj, in such a way as to obtain the linear approximation 

logdet(C-F (51rf2) ^ logdet(Ci + dldO + U (52) 

where we have used that, for X > 0, VxlogdetX = X~^ pp. 641,642]. Then, 
dropping the (irrelevant) constant terms, we end up with the SDP 

minimize tr [(Cj -I- 51^2)”^^] 

subject to C > 0, tr 2 [C] = Id and tri [(p"'" 0 ld)C] = p. (53) 

The minimum of logdet(C -|- 51^2) can be approximated by iteratively solving the 
SDP (53), taking for Cj+i the resulting C of the problem solved with input matrix C*. 
If we set Co = (1 — d)ld 2 as the input for the hrst iteration, the resulting optimization 
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problem coincides with that obtained with the trace heuristic, in which case a rank 
reduction may not occur in the hrst iteration. For successive steps, though, rank 
reduction is generally observed, justifying the initial claim that the log-det heuristic 
is a rehnement of the trace heuristic. 


5.3. Application: Producing X-MEMS with respect to purity 

As an illustration of the effectiveness of the log-det heuristic, let us design economical 
CPTP maps to transform a given input state of purity P into an X-MEMS wrt P, 
denoted by p'p (density matrix obtained by plugging (29), (30) and (31) into (25)). For 


the input state, we take the following X-qubit density matrix of purity P = 1/{N + 1): 

1 


N 


Pn = 


N + 1^ 

k=0 


DtHD 


(64) 


where \D^) are the (totally symmetric) X-qubit Dicke states of k excitations 54,55 
dehned in the computational basis as 


\D^ 


)■■ = 


k\{N-k)\ 


N\ 




k IV—k 


(55) 


with the summation running over every distinct permutation of the sequence of k ones 
and N — k zeros. 


Using SeDuMi, many iterations of the SDP (53) are solved to produce a low rank 


Choi-Jamiolkowski matrix that maps p = pjq top = p'p. Figure]^ shows the evolution 
of rank C* as the iterations progress, where C* denotes the optimal Choi-Jamiolkowski 
matrix found at each step. For completeness and comparison, we also plot the obtained 


optimal values of the SDP (53) at each iteration. Our numerical simulations were run 
with 5 = 0.2 and Co = Id 8) Pp, which corresponds to the CPTP map that collapses 
the entire Bloch ball into the point corresponding to Pp, cf. (48). Plots and 
correspond to the cases N = 3 and = 4, respectively. 

In the case N = 3, plot shows that rank C* = 2 is reached on the 140th iteration. 
From there onwards, every produced C* determines a CPTP map that can be written 
with only two Kraus operator^ff} Remarkably, in this case, the log-det heuristic provides 
minimal Kraus decompositions, since no further decrease can occur (as rank C* = 1 
would correspond to a unitary map). In the case iV = 4, plot|^ shows that rankC* = 3 
is reached after 113 iterations. Although not shown in the plot, we ran 150 iterations 
more and, even so, no further decrease of rankC* was observed. Of course, that does 
not mean that a CPTP map that implement p^ i— p^^ with two Kraus operators does 
not exist, but merely indicates that if it does exist, then it cannot be reached with the 
log-det heuristic. 

We conclude this section by noting that, although the log-det heuristic may not 
always lead to minimal Kraus decompositions, it is at least very effective in producing 


ft See, e.g., 56 Appendix B] for a summary on how to construct sets of Kraus operators for a CPTP 


map from its Choi-Jamiolkowski matrix. 
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a b 




Iteration number 


Iteration number 


Figure 2. Rank decay of the Choi-Jamiolkowski matrix as the iterations of the log- 
det heuristic progress with parameters 6 = 0.2 and Co = Id ® p'p- Plots (a) and (b) 
correspond to the cases TV = 3 and TV = 4, respectively. Tick marks are shown on the 
iterations where a rank drop takes place. For ease of visualization, the optimal values 
of the rank of C and of the objective function of (531 are presented from the third 
iteration onwards. For plot (a), the rank of the initial Choi-Jamiolkowski matrix is 
dr = 8 X 5 = 40, which decays to 24 during the first and second iterations and reaches 
the minimum 2 after 140 iterations. For plot (b), the initial rank is dr = 16 x 8 = 128, 
decaying to 51 during the first and second iterations. After 113 iterations the rank 
reaches 3 and no further decay is observed. 


CPTP maps that implement a desired state transformation with a relatively small 
number of Kraus operators. It should be noted, though, that the efficiency of the 
rank decay scheme is signihcantly limited by the number of qubits involved, as the size 


of problem (53) scales exponentially with N. 


6. Concluding Remarks 

The characterization of multiqubit MEMS is a long-standing open problem in quantum 
information science. In its core, lies the inherent difficulty in quantifying the genuine 
multipartite entanglement of multiqubit systems 57 58 . Notwithstanding, as progresses 


start to be made in this held j9 ,^34 , some preliminar sketches of how iV-qubit MEMS 
look like can be drawn 10 . In this paper, we rely on a recently obtained closed formula 
for the GM-concurrence of A^-qubit X-states to determine — amongst the set of N- 
qubit X-states — those with maximal GM-concurrence for (i) a hxed set of eigenvalues 
and (ii) for a hxed mixedness (as measured by purity), which we refer to as “X-MEMS 
wrt spectrum” and “X-MEMS wrt purity”, respectively. 

Using only elementary algebra, explicit forms of density matrices and maximal 
GM-concurrence were obtained for X-MEMS wrt every possible spectrum. Besides, the 
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unitary transformation that takes an arbitrary iV-qubit state into the corresponding 
X-MEMS wrt spectrum was characterized, generalizing a previous result of Verstraete 


et al. 27 for two-qubits. Although X-MEMS wrt purity had already been identihed 


in we relied on the fact that they form a subset of X-MEMS wrt to spectrum to 
numerically reconstruct them and rigorously prove their optimality via SDP methods. 
Additionally, we formulated as a rank minimization problem the design of quantum 
operations that implements any desired quantum state transformation with the minimal 
number of Kraus operators. Then, applying a heuristic method to specihc examples of 
this optimization problem with N = 3 and 4 qubits, we efficiently characterized low rank 
quantum operations that transform three- and four-qubit states into the corresponding 
X-MEMS wrt purity. 

An extension of our SDP approach to characterize extreme X-states wrt other 
measures of mixedness (e.g. von Neumann entropy) and/or other measures of 
multipartite quantum correlations/nonclassicality (e.g. GM-negativity 30 , global 
quantum discord or the measure introduced in [^) is an interesting line for 
future research. For any desired measures of correlation (c) and mixedness (m), the 
corresponding extreme X-states are, formally, the optimal points of the problem 

maximize c{px) subject to Px > 0, tr px = I, m(px) = nro , (56) 

where mo specihes a desired value of mixedness and the optimization runs over all X- 
density matrices px- Although linearity of c and m would promptly guarantee that 


problem (56) is a SDP, such a form can also be established in certain non-linear cases. 


Indeed, in this paper, we have seen that despite the specihc non-linearities of c (taken 
as the GM-concurrence) and m (taken as the purity), an equivalent SDP was built 
by suitably parametrizing c{px) (cf. 0) and applying a standard trick to turn 
m(px) = nXo into a LMI (cf. Appendix A[ ) . It is thus conceivable that a similar approach 
can handle other choices of non-linear measures c and m. However, determining whether 


(and how) problem (56) can be cast as a SDP is expected to strongly depend on the 
particular choice of measures and, as such, should be considered case by case. 

Along these lines, a particularly interesting problem is to consider a pair of 
optimization problems formed by hxing c as some correlation measure and setting m 
as (i) purity and (ii) von Neumman entropy. The question to be answered here is 
whether the two problems yield the same set of extreme X-states. Such an analysis has 
been conducted by Wei et al. in the bipartite case with c set as concurrence, negativity 
and relative entropy of entanglement [^. Remarkably, different extreme states were 
obtained by changing the mixedness measure, which suggests that the same would occur 
in the multipartite setting. However, an explicit verihcation of this conjecture remains 
an open problem. 
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Appendix A. Converting a quadratic inequality into a linear matrix 
inequality 

In this appendix we briehy review a standard trick to convert quadratic inequality 
constraints into linear matrix inequalities (LMI). In our proof of Theorem]^ this was 
used to write the (convex) quadratically constrained linear program ([M|) in the SDP 


inequality form (35). We illustrate the technique in this particular case. 


First, we use the well-known formula 

^ n \ ^ n n 

+ 2 XkXi , 

^k=l J k=l i>k=l 

to rewrite the quadratic constraint in ([M| as 


(A.l) 


P-l + 2j2^k-2 > 0 (A.2) 

k=i e>k=i 

or, equivalently, 

F-l + 2 j:A-A^Q;'A>0 , (A.3) 

where jn G and G denote the “all-one” n-dimensional vector and matrix, 
respectively, and 

•= (-^15 -^2 , ■ ■ • , An)'*' and := 1„ -|- . (A.4) 

We note that Q~^ is a positive non-singular matrix with eigenvalues 1 ((u — l)-fold 
degenerate) and n -|- 1 (non-degenerate), and its inverse is given by 

1 , 


Qn In 


n + 1 


(A.5) 


Now, we recognize the Ihs of inequality (A.3) as the Schur complement of the 
following (block) matrix of dimension n -|- 1: 


(A.6) 


Qn 

A 

AT 

P-l + 2jJX 


Since Qn > 0, we conclude that (A.3) is equivalent to the constraint that matrix (A.6) 


is positive semidehnite 33, pp. 650-651]. As a result, problem (34) takes the form: 

> 0 (A.7) 




Qn 

A 


maximize — 1 -f 

subject to 

AT 

P — 1 -|- 2jfX 





1 


where b := (2,1,..., I)"*". Up to a constant term and an inversion of sign in the 


objective function (see footnote at page 11), this is equivalent to problem (35). 
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